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Abstract 



A combined fit is made to E791 data onD-» Ktttt, LASS data on 
Kit elastic scattering, and BESII data on J/\& — ► K* (890) Kir. In all 
cases, the Kit S-wave is fitted well with a k pole and A"o(1430); the 
^ | \ k requires an s-dependent width with an Adler zero near threshold. 

\q ; The pole position of the k is at M = (750j^) - i(342 ± 60) MeV. 

CN ■ The E791 collaboration fitted their data using a form factor for the 

production process D — > kit. It is shown that this form factor is not 
needed. The data require pointdike production with an RMS radius 
0^ ■ < 0.38 fm with 95% confidence. 



o 

O . PACS: 13.75Cs, 14.4-OCs, 14.40.Ev 

i/-^ . Keywords: mesons, resonances 

The K in the Kir S-wave has a controversial history and is still a confused 
topic. There are presently definitive data from three processes. 
^ ■ Firstly, LASS data on Kir elastic scattering were fitted originally with an 

effective range formula and no k pole [1]. This fits the data well. 

Secondly, E791 data on D + — > K~n + ii + revealed a peak near threshold 
in the Kir S-wave [2]. A feature requiring the presence of the S-wave is 
^ ■ strong interference with A"* (890); this produces a cos 6 asymmetry in the 

Kit angular distribution. The threshold peak is isotropic and was fitted with 
a conventional Breit-Wigner resonance with mass M = 797 ± 19 ± 43 MeV, 
Y res = 410 ± 43 ± 87 MeV. This corresponds to a pole position M — iF/2 = 
(721 ± 61) - i(291 ± 131) MeV. The data were later re-analysed by Oiler 
using formulae consistent with Chiral Perturbation Theory (ChPT) [3]; he 
reported a k pole at 710— z310 MeV. The CLEO collaboration [4] has reported 
a Dalitz plot in D° — > Kg7t + 1T~ very similar to that of E791. The situation 
is complicated by the appearance of p(770) and other higher p's and also the 
possibility of a — > 7r + 7r~. They were able to reproduce their data with a 
uniform non-resonant Ktc S-wave; this uniform amplitude has a magnitude 
not too far from a broad k resonance. 
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Thirdly, BESII data on J/V -> X*( 890)7^7^ exhibit a low mass Kir 
S-wave peak. It has been fitted with a k pole at 760 ± 20(stat) ± 40(sysi) — 
i(420 ± A5(stat) ± 60(syst)) MeV [5]. A fit by the BES collaboration to the 
same data using a conventional Breit-Wigner parametrisation for the k is 
given in Ref. [6]. This fit gives a k pole at (841 ± 82) - i(309 ± 87) MeV. 

The E791 collaboration has recently refitted the Kit S-wave in their data 
in 37 bins of Kit mass over the entire range from threshold to 1707 MeV 
[7]. In every bin, they determine its magnitude and phase freely, without 
prejudice to how it may be interpreted in terms of k, Kq(1A30) or any other 
component. The objective of the present paper is to show that their data 
are accurately consistent with those from LASS and BES if one introduces 
into the elastic scattering amplitude the Adler zero of ChPT. This combined 
fit eliminates minor inconsistencies between fits to individual sets of data; 
it sharpens the conclusions on the parameters of the k. In passing, it also 
shows that the D decay to kit is point-like within errors. 

The theoretical situation concerning a, k, /o(980) and ao(980) has been 
explored by Oset, Oiler, Pelaez and collaborators [8-11]. They take the real 
parts of S-wave amplitudes from ChPT and obtain the imaginary parts by 
unitarising their amplitudes so as to accomodate rescattering. The key point 
is to explain the difference between elastic scattering and production. Fig. 1 
below shows E791 data, fitted with a variety of form factors for production. 
There is a low mass peak in all cases. BES data shown in Fig. 4(d) also have 
a low mass S-wave peak. However, LASS data on elastic scattering have no 
peak near threshold, Fig. 4(b). Any attempt to fit the production data with 
the LASS amplitude is hopelessly bad. 

Oset et al. show that the reason for the difference is that ChPT demands 
a zero in the elastic scattering amplitude at the Adler point, s = sa = rn 2 K — 
0.5ml ~ 0.23 GeV 2 , not too far from the Ktt threshold at s ~ 0.4 GeV 2 . This 
Adler zero originates from spontaneous Chiral Symmetry breaking, which is 
widely believed to be a key feature of QCD. The k pole lies in the complex 
plane with Re s close to the Kit threshold; the zero distorts the elastic 
scattering amplitude severely and goes a long way towards cancelling the 
nearby k pole. The zero is not present in production processes, where the 
k peak appears undistorted. They report pole positions for the k of 770 — 
i(250 - 425) MeV [9] and 750 - i230 MeV [10]. 

How may these different situations be accomodated into a single formula? 
The assumption which is tested here is that all data may be fitted with an 
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N/D form for the k. The partial wave amplitude for elastic scattering is 
written 

7V(s) 7yexp(2ig)-l 
Tx ^ (s) = ^ = 2* ' (1) 



p(s) is Lorentz invariant phase space 2k/y/s and fc is centre of mass momen- 
tum. The Adler zero at s = sa is included by writing 



9kA s ) = 75 — ^Gexp(- 7 s). (3) 

This is the simplest form which will fit the data. The exponential cut-off 
is a detail required to prevent gj {n increasing indefinitely with s; BES data 
require the amplitude to flatten off or fall at large s. Limitations on the 
parameter 7 will be discussed later. The real part of D(s) is also taken in 
the simplest form which will fit the data, D — s, with D a constant. An 
attempt to fit with D alone fails to fit the observed phase variation with s. 

The denominator D(s) is universal and appears in the amplitude for pro- 
duction. However, N(s) is not universal: it depends on the left-hand cut, 
which is very different for elastic scattering and production. Since left-hand 
singularities for the production process are distant, the amplitude for pro- 
duction data is taken in the form 

_ Aexp(-qg 2 ) 
T D ^ K ^ . (4) 

Here A is a complex constant and the exponential is a form factor depending 
on the momentum q of the k in the D rest frame. 

An important detail is that the Adler zero is a feature of the full Kit S- 
wave. Therefore it is included into the Breit-Wigner amplitudes for i^ (1430) 
and _Ko(1950), which are taken in the form of eqns. (2) and (3). A detail is 
that both resonances are consistent with g 2 Kr> = 0. In the LASS data, there 
are visible resonances at 1430 and 1950 MeV. Both have tails at low mass. To 
keep the amplitude unitary below the inmelastic threshold, the Dalitz-Tuan 
prescription [12] is used. The total S-matrix is written as the product of Si 
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Ktt mass (GeV) Krr mass (GeV) 

Figure 1: Fits to the magnitude of the k amplitude in E791 data from Ref. 
[7] for four values of a in the form factor; a is in units of GeV -2 . 



for the three individual resonances: 

Stot = S1S2S3, (5) 
T = (S-l)/(2i), (6) 

For decays of D and J/^, there are large numbers of open channels, none of 
which is close to the unitarity limit. In this case, unitarity plays no role and 
the standard isobar-model is adopted, adding amplitudes from eqn. (4) with 
fitted complex coupling constants A, for each resonance. 

In Ref. [7], E791 adopted the form factor F = exp(— aq 2 ) for the produc- 
tion process D — > (Ktt)tt with a = 2.08 GeV -2 . Here this form factor will 
be adjusted to obtain a fit consistent with all sets of data. Fig. 1 shows the 
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fit to E791 data for four values of a. Panel (a) gives the best fit with a = 0. 
In (b) - (d), the exponent of the form factor is increased in equal steps to 
the E791 value in (d). Fig. 2(a) shows x 2 as a function of a. There is an 
optimum just below zero for both magnitude and phase, but consistent with 
a = within the errors. The RMS radius corresponding to this form factor 
is zero and < 0.38 fm with 95% confidence. 



x 




2 2 

a (GeV~ 2 ) a (GeV~ 2 ) 



Figure 2: \ 2 f° r E791 data as a function of a in the exponential form factor: 
(a) for the magnitude of the amplitude, (b) for the phase. 

Fig. 3 shows corresponding fits to the phase of the k. The sensitivity to 
a is less evident to the eye; nonetheless it again optimises at zero within the 
errors, as shown on Fig. 2(b). On careful inspection of Fig. 3(d), the fit is 
systematically high at both ends of the mass range. 

A bin-by-bin fit to BESII data is reported in Ref. [5]. Magnitudes and 
phases from that analysis, shown in Figs. 4(c) and (d), are fitted here si- 
multaneously with LASS and E791 data. Variations of k magnitude and 
phase with Kit mass are consistent between LASS and BES II data within 
the errors. 

An important ingredient is the mass and width of K (1A30). On Figs. 1 
and 3, it is obviously needed by E791 data, but a fit to these data alone shows 
considerable flexibility in its parameters. The same is true for LASS data. 
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Figure 3: Fits to the phase of the k amplitude in E791 data for the same a 
as Fig. 1. 



In BES data, the K (1430) in the channel K* (890) K ( 1430) is a large signal 
with well defined centroid and width. Those data constrain the mass and 
width of this resonance. A complication is that there is a smaller X 2 (1430) 
signal at the same place, but it is well separated by the amplitude analysis. 
Fig. 4(c) shows the mass projection of the 1430 MeV Kn peak in BES data, 
allowing an intensity ratio K 2 (U30) / K (U30) = 0.75 [5]. 

Now we come to numerical details of the fit. The Kq(1A30) amplitude 
is sensitive to its coupling to Kr]'. Above this threshold, the phase of the 
amplitude varies rapidly with mass. The opening of the Krj' channel also 
affects the line-shape. It turns out that E791 data have only weak sensitivity 
to the ratio r = g 2 (Krj') / g 2 (Kir). If fitted alone, they optimise with zero 
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Kn mass (GeV) Kn mass (GeV) 



Figure 4: Fits to (a) the phase, (b) the magnitude of LASS amplitudes for 
elastic scattering, (c) the phase, (d) the magnitude of BES data, (e) the 1430 
MeV peak in the Kir mass projection of BES data, (f) the phase in E791 
data. 
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coupling to Ki] f , but x 2 changes only by 3.6 for a ratio r — 1. LASS data 
alone optimise at r = 1.30 ± 0.16. The full partial wave analysis of BES 
data is quite sensitive to r through interferences with other partial waves 
and optimises at r = 1.0 ± 0.2. We adopt the mean value r = 1.15, with an 
estimated error of ±0.22. The sensitivity of different sets of data to different 
components highlights the benefits of fitting several sets of data together. 

Some constraint is needed on the parameter 7 of the exponential in eqn. 
(4). It tends to run away to a large value. However, this leads to a rapid 
reduction in the width of the k, at the top end of the mass range and to an 
increase in the k intensity above 1.6 GeV. That range is not covered by the 
bin-by-bin data being fitted here, but is available in the full partial wave 
analysis. Those data set upper limits on the magnitude of the k component 
between Kir masses of 1.6 and 2.2 GeV and rule out any significant increase 
in the n amplitude in this range. We therefore constrain the parameter a to 
be < 0.25 to avoid such an increase. 

The fit to LASS data is not perfect in the mass range 1.15 to 1.35 GeV, 
see Figs. 4(a) and (b). It is possible that this is due to some decay to Knun, 
presently not known experimentally; the decay to k<j is a possibility which 
appears naturally in the linear o model. Some improvement to the fit is 
possible by replacing G of eqn. (4) by G(l + (3s). However, /3 then goes 
negative and again creates a problem in the fit to BES data in the mass 
range 1.6 to 2.2 GeV. The negative value of (3 produces a zero either in this 
mass range or close to it and a large increase in the n amplitude. Again this 
is ruled out by the BES data. So we set (3 = 0. 

The mass and width of -Ko(1950) are varied in the range of values quoted 
for solution A by LASS. This solution is favoured by the analysis of Biittiker, 
Descotes-Genon and Moussallam [13]. The phases of the highest three / = 
1/2 points of LASS (1.53 GeV upwards) cannot be fitted accurately with any 
variation of parameters and are omitted from the final fit; these points have 
little sensitivity to g 2 Kri , or parameters of fTo(1950). For the coherent sum 
of / = 1/2 and 1 = 3/2 amplitudes, LASS found results lying outside the 
unitarity circle, so it is possible that there is a significant sensitivity to the 
/ = 3/2 component which was subtracted by LASS. 

Equally good fits may be obtained to all data sets for a large range of D 
in eqns. (2) and (3). The only real requirement is that y/Do is well above 
the mass of i^o(1430), since the phase of the k rises slowly to ~ 72° at 1.7 
GeV. The value D = 3.3 2 GeV 2 is adopted from earlier fits to BES data. 
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The optimum fit gives G = 24.5 GeV 2 and g 2 (Kr])/ g 2 (Kn) = 0.06 ± 0.02. 
The fit is insensitive to ry = g 2 (Ki]')/ g 2 (Kii) for the k because of cross-talk 
with _K" (1430) — > Kvj '. A range of values for r v > is possible from to 0.57 
with little change in x 2 - The pole position of the k is rather sensitive to this 
variation. It varies from 743 — i367 MeV to 757 — i317 MeV. We adopt a 
mean value (750^5°) — i(342 ± 60) MeV. Errors are mostly systematic and 
depend on variations of the parameter a of the exponential cut-off for the k 
and r = g 2 (Kri')/g 2 (Kir) for K (U30) and k. 

For X o (1430), fitted parameters are M = 1.517 GeV, g 2 (Kn) = 0.353 
GeV 2 for the chosen value g 2 (Ki]')/g 2 (Ki]) = 1.15. The pole position is at 
M = (1406 ±29) -i(175±20) MeV. Data on Kir -> Krj' would reduce these 
systematic errors. There is the strong possibility that either or both of k and 
-£To(1430) decay to Kirnir, so data on that channel would be helpful. 

Some comment is required on the phase of Kn elastic scattering. How 
is it possible to have a pole without the phase going through 90° in the 
mass range from threshold to 1.4 GeV? The clue is that the pole lies at 
s = 0.445 — i0.513 GeV 2 ; its real part lies almost below the Ktt threshold at 
s = 0.401 GeV 2 . A Breit-Wigner resonance of constant width would have a 
phase of 90° at a mass of 667 MeV. The effect of p(s) is to produce a phase 
rising rapidly from zero at threshold. It is essential to realise that the phase 
varies rapidly as one goes off the real s-axis. There is a rapid phase variation 
from the imaginary part of (s — sa) and also due to p(s) becoming complex. 
Near the pole, the phase is 90° away from its threshold value on the real 
s-axis; this was pointed out by Oiler in Ref . [3] . The combined effect of the 
Adler zero and p(s) is to retard the phase on the real s-axis by roughly 90° 
from the value produced by the pole alone. This is a major distortion from 
the familiar resonance of constant width. 

Although this may appear surprising, it succeeds in fitting LASS, E791 
and BES data consistently in magnitude and phase. It is also consistent 
with ChPT. Adding one more power of s to either or both of the real and 
imaginary parts of D(s) in eqn. (2) does not remove the pole; the pole 
position changes within the quoted errors. The reason for this stability is 
that LASS and E791 data plus CHPT determine the phase all the way from 
the Adler zero to ~ 1500 MeV. Zheng et al. [14] have examined fits to LASS 
data using a range of formulae departing from strict ChPT. Their conclusion 
is that a pole is required if the scattering length is ao < 0.34m" 1 . From the 
fit reported here, ao = 0.23 ± 0.04m" 1 , compared with the value 0.19m 7i : 1 of 
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ChPT at order p 4 [15]. 

In conclusion, a combined fit to E791, LASS and BES data separates the 
k and X (1430) components of the Kn S-wave and provides a consistent fit to 
all three sets of data. They all display the same phase variation as a function 
of mass for the k component. No 'background' is required within present 
errors. If such background is present, it is the same in all three sets of data 
and can be absorbed algebraically into the s-dependence of the amplitude. 
The form factor for production of the Kir S-wave in E791 data optimises 
at F — 1, corresponding to point-like production. The RMS radius for the 
production process is < 0.38 fm with 95% confidence. 

I wish to thank Prof. Brian Meadows for careful and helpful suggestions. 

References 

[1] D. Aston et al. (LASS), Nucl. Phys. B296 (1988) 493. 

[2] E.M. Aitala et al. (E791), Phys. Rev. Lett. 89 (2002) 121801. 

[3] J.A. Oiler, Phys. Rev. D71 (2005) 054030. 

[4] S. Kopp et al., (CLEO) Phys. Rev. D63 (2001) 092001. 

[5] D.V. Bugg, Euro. Phys. J A25 (2005) 107. 

[6] M. Ablikim et al. (BES II), hep-ex/0506055. 

[7] E.M. Aitala et al. (E791), hep-ex/0507099, submitted to Phys. Rev. D. 

[8] J.A. Oiler and E. Oset, Nucl. Phys. A620 (1997) 438. 

[9] J.A. Oiler, E. Oset and J.R. Pelaez, Phys. Rev. D59 (1999) 074001. 

[10] J.A. Oiler and E. Oset, Phys. Rev. D60 (1999) 074023. 



10 



[11] A. Gomez Nicola and J.R. Pelaez, Phys. Rev. D65 (2002) 054009. 

[12] R.H. Dalitz and S. Tuan, Ann. Phys. (N.Y.) 10 (1960) 307. 

[13] P. Buttiker, S. Descotes-Genon and B. Moussallam, Euro Phys. J C33 
(2004) 409. 

[14] H.Q. Zheng et al, Nucl. Phys. A733 (2004) 235. 

[15] V. Bernard, N. Kaiser and U.G. Meissner, Nucl. Phys. B357 (1991) 129. 



11 



